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We incorporate the effective restoration of (7(1) a symmetry in the 2+1 flavor entanglement Polyakov-loop 
extended Nambu-Jona-Lasinio (EPNJL) model by introducing a temperature-dependent strength K{T) to the 
Kohayashi-Maskawa-’t Hooft (KMT) determinant interaction. T dependence of K{T) is well determined from 
pion and ao-meson screening masses obtained hy lattice QCD (LQCD) simulations with improved p4 stag¬ 
gered fermions. The strength is strongly suppressed in the vicinity of the pseudocritical temperature of chiral 
transition. The EPNJL model with the K{T) well reproduces meson susceptibilities calculated by LQCD with 
domain-wall fermions. The model shows that the chiral transition is second order at the “light-quark chiral- 
limit” point where the light quark mass is zero and the strange quark mass is fixed at the physical value. This 
indicates that there exists a tricritical point. Hence the location is estimated. 

PACS numbers: 11.30.Rd, 12.40.-y, 21.65.Qr, 25.75.Nq 


I. INTRODUCTION 


Meson masses are important quantities to understand the 
properties of Quantum Chromodynamics (QCD) vacuum. 
For example, the difference between pion and sigma-meson 
masses is mainly originated in the spontaneous breaking of 
chiral symmetry, so that the restoration can be determined 
from temperature (T) dependence of the mass difference. 
Similar analysis is possible for the effective restoration of 
[7(1)A symmetry through the difference between pion and ao- 
meson (Lorentz scalar and isovector meson) masses. 

(7(1) A symmetry is explicitly broken by the axial anomaly 
and the current quark mass. In the effective model, the U (1)a 
anomaly is simulated by the Kob^ashi-Maskawa-’t Hooft 
(KMT) determinant interaction |lTl|2]- The coupling constant 
K of the KMT interaction is proportional to the instanton den¬ 
sity screened by the medium with finite T Q. Hence K be¬ 
comes small as T increases: K = K{T). Pisarski and Yaffe 
discussed the suppression S{T) = K(T)/K{0) for high T, 
T ^ ^Tc for the pseudocritical temperature Tc of chiral 
transition, by calculating the Debye-type screening lU : 


S{T) = exp -Tt^p^T^[ 


:N, 


:Ni 


( 1 ) 


where Nc (Nf) is the number of colors (flavors) and the typ¬ 
ical instanton radius p is about 1/3 fm, and hence the sup¬ 
pression parameter b is about O.TOTc for Nc = Nf = 3 of 
our interest i] ; note that 2+1 flavor LQCD simulations show 
Tc = 154 ± 9 MeV lUl-Ql- This phenomenon is called “effec¬ 
tive restoration of U{1)a symmetry”, since (7(1)a symmetry 
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is always broken in the current-operator level but effectively 
restored at higher T in the vacuum expectation value. 

Figure[T]shows the current status of knowledge on 2+1 fla¬ 
vor phase diagram for various values of light-quark mass m; 
and strange-quark mass mg. QCD shows a first-order phase 
transition associated with the breaking of chiral (Zf) symme¬ 
try at the lower left (upper right) corner lH, . When mi and 
TOs are finite, these first-order transitions become second or¬ 
der of 3d Ising {Z{2)) universality class, as shown by the solid 
lines |[1,|1]. 

However, the order of chiral transition and its universal¬ 
ity class are unknown on the vertical line of mi =0 and 
TOs > 0, and it is considered to be related to the effec¬ 
tive ?7(1)a restoration. In the two-flavor chiral limit of 
{mi, ms) = (0,oo) at the upper left corner, for example, 
the order may be second order of 0(4) class if the effective 
restoration is not completed at T = Tc, because the chiral 
symmetry becomes SUi^{2) x S'(7r(2) isomorphic to 0(4) in 
the situation and the transition is then expected to be in the 3d 
0(4) universality class |[1,|^. When (7(1)a symmetry is re¬ 
stored completely at T — Tc, it was suggested in Ref. that 
the chiral transition becomes the first order. Recently, how¬ 
ever, it was pointed out in Ref. iH that the second order is 
still possible. In this case, the universality class is not 0(4) 
but (7l(2) X (7r,(2). There are many lattice QCD (LQCD) 
simulations made so far to clarify the order and its universal¬ 
ity class in the two-flavor chiral limit of (to/, TOs) = (0, oo) 
and the light-quark chiral limit where to/ vanishes with TOs 
fixed at t he p hysical value, but these are still controversial; 
see Refs. and therein. 

Very recently, the effective restoration of U{1)a symmetry 
was investigated by pion and ao-meson screening masses cal¬ 
culated with LQCD with improved p4 staggered fermions 
and also by meson susceptibilities calculated with LQCD with 
domain-wall fermions 1251] . The effective restoration of 
[/(1)a symmetry thus becomes an important current issue. 

In LQCD, pole and screening masses are evaluated from 
the exponential decay of mesonic correlation functions in the 
temporal and spatial directions, respectively, but for finite T 
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OO 2 flavor pure gauge 



meson screening mass scr is difficult, where ^ denotes the 
typ e of meson. In fact, only a few trials were made so far ifs^ 
1^. The first problem is the regularization. The regularization 
widely used is the three-dimensional momentum cutoff, but it 
breaks Lorentz invariance at T = 0 and spatial-translation 
invariance at any T. This generates un phy sical oscillations in 
the spatial correlation function 7744 (r) js!]]. This refuses us to 
determine scr from the asymptotic form of 7744 (r) as 


= - to (2) 

r^oo dr 

This problem can be solved llsll] by introducing the Lorentz- 
invariant Pauli-Villars (PV) regularization lls^ . 


Fig. 1: A schematic phase diagram of 2-1-1 flavor QCD as a func¬ 
tion of light-quark mass mi and strange-quark mass ms . A tricritical 
point is likely to appear on the ms axis; the location is shown by 
{mi, ms) = (0, mr^). The solid lines stand for second-order tran¬ 
sitions belonging to the universality class labeled, where the labels 
Z{2) and 0(4) mean the 3d Ising and the 3d 0(4) class, respec¬ 
tively. 


Even after the unphysical oscillations are removed, heavy 
numerical calculations are still required to obtain r]^^{r) at 
large r il. This is the second problem. In the model calcula¬ 
tion, the spatial correlation function is obtained first in the mo¬ 
mentum representation {q = ±|q|) as Hence we 

have to make the Fourier transform from (0, q^) to 77 ^^ (r): 

dq qx««(0, (3) 


the lattice size is smaller in the temporal direction than in the 
spatial direction. This makes LQCD simulations less feasi¬ 
ble for pole masses than for screening masses. The problem 
is getting serious as T increases. This is the reason why me¬ 
son screening masses are calculated in most of LQCD sim¬ 
ulations. In fact, as mentioned above, state-of-the-art LQCD 
calculations were done for meson screening masses with large 
volumes (16^ x 4, 24^ x 6 , 32^ x 8 ) in a wide range of T 1^ . 

The effective model is suitable for qualitative understand¬ 
ing of QCD. In fact, the QCD phase diagram, the properties 
of light-meson pole masses and the T dependence of U{1)a 
anomaly have been studied extensively with the Poly^ov- 
loop extended Nambu-Jona-Lasinio (PNJL) model 
The PNJL model can treat the deconfinement and chiral tran¬ 
sitions simultaneously, but cannot reproduce their coincidence 
seen in 2-flavor LQCD data, when the model parameters are 
properly set to reproduce Tc calculated with LQCD ll^ . This 
problem is solved by introducing the four-quark vertex de¬ 
pending on the Polyakov-loop |4^ . The model with the 
entanglement vertex is refer to as the entanglement-PNJL (EP- 
NJL) model. The EPNJL model is successful in reproducing 
the phase diagram in 2-flavor QCD at imaginary chemical po¬ 
tential |4N] and real isospin chemical potential and 
well accounts for the phase diagram in the mi-rris plane, i.e., 
the so-called Columbia plot in 2 h- 1 flavor QCD 1^ . So far, 
T dependence of t/(l)A anomaly and low-lying meson pole 
masses was studied extensively by both the PNJL and EPNJL 
models ll^l^lddll^ . but T dependence of meson screen¬ 
ing masses was investigated by the EPNJL model only in our 
previous work . 

In NJL-type models, it is well known that the calculation of 


The q integration is quite hard particularly at large r, since the 
integrand consists of a slowly damping function 
and a highly oscillating function If XJ5(0! 9^) has a pole 
below the cut in the complex q plane, one can easily deter¬ 
mine Mf snr from the pole location. In the old formulation 
of Ref. iDin . the condition was not satisfied, since logarith¬ 
mic cuts appear in the vicinity of the real q axis in addition 
to physical cuts. Very recently we solved the problem in our 
previous paper S, showing that the logarithmic cuts near the 
real q axis are unphysical and removable. In the new formu¬ 
lation based on the PV regularization, there is no logarithmic 
cut and a pole appears below physical cuts, as shown later. 

In this paper, we incorporate the effective restoration of 
f7(l)A symmetry in the 2 h- 1 flavor EPNJL model by intro¬ 
ducing a T-dependent coupling strength K{T) to the KMT 
interaction. T dependence of K{T) is well determined from 
state-of-the-art 2 h- 1 flavor LQCD results ll2^ on pion and ag- 
meson screening masses. For the derivation of meson screen¬ 
ing mass, we extend the previous prescription of Ref. 
for 2 flavors to 2 h- 1 flavors. The K{T) determined from the 
LQCD data is strongly suppressed near Tc- Using the param¬ 
eter set, we show that the chiral transition is second order in 
the light-quark chiral limit. This result indicates that there ex¬ 
ists a tricritical point near the “light-quark chiral-limit” point 
in the mi-mg plane. We then estimate the location. 

We recapitulate the EPNJL model and the method of calcu¬ 
lating meson screening masses in Sec. |II]and show the results 
of numerical calculations in Sec. |III] Section|IV]is devoted to 
a summary. 
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II. MODEL SETTING 


A. EPNJL model 


We start with the 2+1 flavor EPNJL model l43l] . The 
Lagrangian density is 


£ - mo)V' + Gs{^) '^[{tpXaip)'^ + 


a—0 


- K(T) det'!/'/(l+75)'*/'/' + det' 0/(1 - 75)1/'/' 
-Ui<P[A],^A],T) 


(4) 


with quark fields 0 = ( 0 u, 0 d, 0 s)^ and D'^ = d'^ + iA’' 
with A’' = S^g(A°)ata /2 = - 6 ^ig{A 4 )ata /2 for the gauge 
coupling g, where the Xa (fa) are the Gell-Mann matrices in 
flavor (color) space and Aq = a/2/3 I for the unit matrix I in 
flavor space. The determinant in (|4|l is taken in flavor space. 
For the 2+1 flavor system, the current quark masses mg = 
diag(TO„, Ws) satisfy a relation > m; = = m^. 

In the EPNJL model, the coupling strength Gs{‘P) of the 
scalar-type four-quark interaction depends on the Polyakov 
loop <I> and its Hermitian conjugate 0 as 


with L = exp[*A 4 /T] = exp[idiag(A/, for 

real variables A/ satisfying A^ + A^^ + A^^ = 0. For zero 
quark chemical potential where = ^, one can set ^ 4 ^ = 0 
and determine the others as A^ = —A^ = cos“^[(3t? — 

l)/2]. 

We use the logarithm-type Polyakov-loop potential of 
Ref. ll^ as U. The parameter set in U has already been de¬ 
termined from LQCD data at finite T in the pure gauge limit. 
The potential has a parameter Tq and yields a first-order de¬ 
confinement phase transition at T = Tg. The parameter used 
to be set to Tg = 270 MeV, since LQCD data show the phase 
transition at T = 270 MeV in the pure gauge limit. In full 
QCD with dynamical quarks, however, the PNJL model with 
this value of Tg is found not to explain LQCD results. Nowa¬ 
days, Tg is then rescaled to reproduce the LQCD results. In 
the present case, we take Tg = 180 MeV so that the EP¬ 
NJL model can reproduce LQCD results for the pseudocriti- 
cal temperature T/econf deconfinement transition; actually, 
ydeconf ^ jjj jgg EPNJL model and 170 ± 7 MeV 

in LQCD lUl. 

Making the mean field approximation (MFA) to dUi leads to 
the linearized Lagrangian density 

C^^^ = i^S-^i^-UM-U{<P[A],§[A],T) ( 8 ) 


Gs(<J>) = Gs(0) [1 - ai<?0 - a2(<?^ + <?^)] . (5) 

This entanglement coupling is charge-conjugation and 
symmetric. When ai = a 2 = 0, the EPNJL model is reduced 
to the PNJL model. We set a 2 = 0 for simplicity, since the a 2 
term yields the same effect as the ai term in the present anal¬ 
ysis. As shown later in Sec. |III] the value of ai is determined 
from LQCD data on pion and og-meson screening masses; the 
resulting value is ai = 1 . 0 . 

For T dependence of K (T), we assume the following form 
phenomenologically: 


K{T) = 


K{Q) (T < Ti) 

(T > Ti) 


( 6 ) 


For high T satisfying T ^ Ti, the form (| 6 ]l is reduced to O- 
As shown later in Sec. mil the values of Ti and b are well 
determined from LQCD data on pion and og-meson screening 
masses; the resulting values are Ti = 0.79Tc = 121 MeV and 
b = 0.23Tc = 36 MeV. 

After the Pisarski-Yaffe discussion on S{T), T dependence 
of the instanton density was estimated theoretically by the the 
instanton-liquid model |@ 1 , but the estimation is applicable 
only for T > 2Tc. For this reason, in Ref. 14^ . a Woods- 
Saxon form (1 + j-i with two parameters T[ and b' 

was used phenomenologically for K{T)/K{0). The present 
form (| 6 ]l has T dependence similar to the Woods-Saxon form. 

In the EPNJL model, the time component of A^ is treated 
as a homogeneous and static background field, which is gov¬ 
erned by the Polyakov-loop potential U. In the Polyakov 
gauge, and 3 are obtained by 


with the quark propagator 

,5 = (n,0"-7gAO-M)-/ (9) 

where M = diag(M„, Mj,, Mg) with 

Mu = mu - 4Gs(T)(t„ -I- 2K{T)(Jd(Js, 

Md = md- 'iGs{<P)ad + 2K{T)asau, 

Ms = ms - 4Gs(T)(Ts -I- 2K{T)au<7d, 

and ( 7 f means the chiral condensate for flavor /. The 

mesonic potential Um is 

Um = 2 Gs{^){al + al + al) - 4iT(T)a„ada«. 


Making the path integral over quark fields, one can get the 
thermodynamic potential (per unit volume) as 


17 = [/m+G-2 

f—u,d^s * 


d^v r 

(27r)3 : 


3T, 


'pj 


-I-— In [1-f 3(T-f + e 


-f - In [ 1 -f 3(t?-f ^^P’^)e ^^P'f + e 


( 10 ) 


with Epj = + M'j and (3 = 1/T. We determine the 

mean-field variables (X = ai,as, 3) from the stationary 
conditions: 



( 11 ) 


<J>=itr/T), 0 = itr/T*) 


(7) 


where isospin symmetry is assumed for the light-quark sector, 
i.e., (Ji = ( 7 u = Cd. 
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On the right-hand side of (fTOl) . the first term (vacuum term) 
in the momentum integral diverges. We then use the PV reg¬ 
ularization llsills^l . In the scheme, the integral I{Mf,q) is 
regularized as 


2 

= ( 12 ) 

a—0 

where Mf.Q = Mf and the Mf.^a. (ct > 1 ) mean masses of 
auxiliary particles. The parameters and Ca should sat¬ 
isfy the condition = 0 - We 

then assume (Co,Ci,C 2 ) = ( 1 , 1 ,— 2 ) and Mj.^) = 

{Mj + 2yl^, Mj + yl^). We keep the parameter yl finite even 
after the subtraction (fTZt . since the present model is non- 
renormalizable. The parameters are taken from Ref. ll^ and 
they are mi = 6.2 MeV, m, = 175.0 MeV, Gs(0)yl2 = 2.35 
and K(Qi)A^ = 27.8 for A = 795 MeV. This parameter set 
reproduces mesonic observables at vacuum, i.e., the pion and 
kaon decay constants (/tt = 92 MeV and Jk = 105 MeV) 
and their masses (Mjr = 141 MeV and Mk = 512 MeV) 
and the ry'-meson mass (M^/ = 920 MeV). In the present 
work, we analyze LQCD results of Ref. for pion and oq- 
meson screening masses. In the LQCD simulation, the pion 
mass at vacuum (T = 0) is 175 MeV and a bit heav¬ 

ier than the experimental value 138 MeV. We then change 
mi to 9.9 MeV in the EPNJL model in order to reproduce 
M,n.(0) = 175 MeV. This parameter set yields Ma(,(0) = 711 
MeV, Mr^iQ) = 481 MeV and M^{Q) = 537 MeV as ao,?? 
and cr meson pole masses at vacuum. 


B. Meson pole mass 

We derive the equations for pion and oo-meson pole 
masses, following Ref ll^ 1^ . The current corresponding 
to a meson of type ^ is 

Jj(x) = '0(x)/5V’(x) - (13) 

where /V = for tt meson and Fa^ = A 3 for OQ-meson. 

We denote the Fourier transform of the mesonic correlation 
function 77 jj(x) = (0|T J^(x) j|(0)) |0) by as 

= Xid^o^f) = * J d‘‘xe*«'^ry^^(x), (14) 

where q = ±|q| for q = (qg, q) and T stands for the time- 
ordered product. Using the random-phase (ring) approxima¬ 
tion, one can obtain the Schwinger-Dyson equation 

Xii = -^44 + 2 ^ (15) 


with p' = (po + iA 4 ,p), where the trace Tr is taken in 
flavor, Dirac and color spaces. Here the quark propagator 
S{p) in momentum space is diagonal in flavor space: S{p) = 
diag{Su,Sd,Ss)- For ^ = tt and og, furthermore, and 
iTjj/ are diagonal (G^^/ = because 

we impose isospin symmetry for the light-quark sector and 
employ the random-phase approximation. One can then eas¬ 
ily solve the Schwinger-Dyson equation for ^ = tt and og: 


= 

1 - 2 G^n^ 


with the effective couplings G^n- and Gao defined by 
Ga„ = G.(<?) + 

G^ = G,(<?) - \K{T)as. 


(17) 


(18) 

(19) 


As for r = 0, 77^ and TTaj, have the following explicit 
forms: 

nao = ^y^(A3)/'/(A3)//' 
fj' 


dAp 


trc 


{idp'+ + + Mf,) 


J (2^)4 

= 4i[h+l2-{q"^-AMdh] 

Utt = ^y^(A3)/'/(A3)//' 


{{p' + q)2 - M2}(p'2 _ M],) 


fJ' 


dAp 


trc 


(* 75 ) 


{iii{p' + (lY + ^f} 


(2^)4 c.a^v 

{h^p''" + ixlf,)- 


x(*75) 


= 4z[/i +I 2 - q^h], 


and 


h = 
h = 
h = 


trc 


d^p 

(27r)4 “^np'2 - mV 
d'^p 


trc 


(27r)4 L(p'-I-q)2 — m2. 
d'^p 


trc 


( 20 ) 


(27r)4 ‘^L{(p'-fq)2-M2}(p'2-M2)J’ 


( 21 ) 

( 22 ) 

(23) 

(24) 


where trc,d (trc) means the trace in color and Dirac spaces 
(color space) and M = Mu = Md- For finite T, the corre¬ 
sponding equations are obtained by the replacement 


pg —>■ iuJn = i{2n + l) 7 rT, 



for where G^/^" is an effective four-quark interaction and 
is the one-loop polarization function defined by 

7T^e(q2) = (-Q J ^Tr{r^tSip' + q)r^nS{p'))il6) 


The meson pole mass Mj^poie is a pole of XiiYYd^) 
the complex qg plane. Taking the rest frame q = (qg, 0) for 
convenience, one can get the equation for poie as 


[l-2G57T5(qg^0)]l 


QO —A^^ipole ^"2 


(26) 
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where F is the decay width to qq continuum. The method of 
calculating meson pole masses is well established in the PNJL 


model pfilL55D. 


C. Meson screening mass 


We derive the equations for pion and ag-meson screen¬ 
ing masses, following Ref. This is an extension of the 
method of Ref. (4^ for 2 flavors to 2 h- 1 flavors. 

As mentioned in Sec. H] it is not easy to make the Fourier 
transform from ^ particularly at large r. 

When the direct integration on the real q axis is difficult, one 
can consider a contour integral in the complex q plane by us¬ 
ing the Cauchy’s integral theorem. However, X€€(0^^) has 
logarithmic cuts in the vicinity of the real q axis 11 ^ , and 
it is reported in Ref. m that heavy numerical calculations 
are necessary for evaluating the cut effects. In our previous 
work ll4^ . we showed that these logarithmic cuts are unphys¬ 
ical and removable. Actually, we have no logarithmic cut, 
when analytic continuation is made for the I 3 (q) after p inte¬ 
gration. Namely, the Matsubara summation over n should be 
taken after the p integration in (l25T l. We then express Jg®® as 
an infinite series of analytic functions: 


Nc oa 2 

I^{o,e) = ^TY^ E 

j — l n— — oo Q;=0 

f d^p r 1 1 


iT 

27r^ 

iT 

Fnq 


(27r)3 Lp2 (p + q)2 




dx / dk- 


j,n,a. 


Ca sin ^ 


0 Jo + ix- x^)q^ + 

(27) 


j,n^a 


4 ' j,n,a 


with 


M,, „,„(T) = Y/M2 + {(2n+l)7rT + Af}2, (28) 

where Ma = Mu;a = Md-,a- We have numerically con¬ 
firmed that the convergence of n-summation is quite fast in 
(|27]l. Each term of l^‘^^{0,q^) has two physical cuts on 
the imaginary axis, one is an upward vertical line starting 
from q = 2 iMj^n,a and the other is a downward verti¬ 
cal line from q = —2iMj^n,a- The lowest branch point is 
q = 2iMj^i n=o,a=o- The value is the meson screening mass 
in which there is no interaction between a quark and an anti¬ 
quark, i.e., = 0. Hence we may call 2Mj=i „=o,a=o “the 

threshold mass”. 

We can obtain the meson screening mass Mj^scr as a pole 

[1 - 2G£77.(0, f __=0. (29) 

If the pole at q = is well isolated from the cut, i.e., 

T7{,scr < 2,Mj=i^n=o,a=07 one Can determine the screening 
mass from the pole location without making the q integral. In 
the high-T limit, the condition tends to scr < 27rT. 


D. Meson susceptibility 


We consider meson susceptibilities for ^ = 7r,ao,p 
and a. In LQCD simulations of Refs. iSlUt], the are 
defined in Euclidean spacetime xe = (t, x) as 

xr = IJ {JdF x)4 (0, 0)). (30) 


In the simulations, J^- and J^i are assumed to 
have no s-quark component for simplicity: namely, 

= _J2f=u,dj’f'J’f ~ i'l2f=u,dj’f'J’f) and = 

E f=u,d i’Flodf - (E f=u,d V’/*75V'/)- For consistency, we 
take the same assumption also in the present analysis, and 
denote the mesons with no s-quark (ns) component by CThs 
and pns- The factor 1/2 is introduced to define the XSUS 
single-flavor quantities. 

The x|“® is related to the the Matsubara Green’s function 
xf^( 94 : in the momentum representation as 

xr = ^x|(7l,q")|,,=o,q=o> (31) 


and x^^ is obtainable from (fTTl l for tt and oq mesons. Eor 
Pns meson, we have to consider a mixing between and 
Vs = dsdods- As a result, one can obtain ns 11551] 


Xr;„ 


(1 


)n. 


det [/ - 2G7T] 
where I is the unit matrix and 


G = 

for the elements 


Gri^ri^ Gri^nr,^ 


Gr]naT]B 


77 = 


77„ 


0 77, 


77 ns T^ns 


(32) 


(33) 


Gr^BVs = Gs(<7>), (34) 

= Gs(<7) + lK{T)as, (35) 

Gr/aJ7ns = = —^K{T)ai. (36) 


In the isospin symmetric case we consider, the polarization 
functions and have the same function form as 


= n^riMs), (37) 

Gkrinarins — ^[-^(141), (38) 

where note that 7T,r (-Tfs) is a function of not M but Mg. Sim¬ 
ilarly, Xo-naCTna is Obtainable from ( l32l i with K{T) replaced by 
-K{T) and 77^ by 


III. NUMERICAL RESULTS 
A. Meson screening masses 

The EPNJL model has three adjustable parameters, ai in 
the entanglement coupling Gs (^) and b and Ti in the KMT 
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Fig. 2: T dependence of pion and ao-meson screening masses, 
M^,scr and MaQ,aci- The solid (dot-dash) line denotes 
(Mag,aci) calculated by the EPNJL model, whereas the dotted line 
corresponds to the threshold mass. LQCD data are taken from 
Ref. 12311 : closed squares (open circles) correspond to the 2+1 fla¬ 
vor data for M.,r,scr (Mag, bci)- In Ref. 1^ . Tc was considered 
196 MeV, but it was refined to 154 ± 9 MeV The latest 

value is taken in this figure. 


interaction K{T). The se p arameters can be clearly deter¬ 
mined from LQCD data 11^ for pion and oo-meson screening 
masses, ^cr and Mag^sa, as shown below. 

Figure |2] shows T dependence of Mtt, scr and Mag^sa- Best 
htting is obtained, when ai = 1.0, Ti = 0.79Tc = i21 MeV 
and b = 0.23Tc = 36 MeV. Actually, the EPNJL results 
(solid and dot-dash lines) with this parameter set well ac¬ 
count for LQCD data for both MTr^scr and scr- The 
parameters thus obtained indicate the strong suppression of 
K(T) in the vicinity of Tg. The mass difference AMsci{T) = 
Mag^scr{T) — M^_scr(T) is Sensitive to K(T) because of (fTSl l 
and ( fT9] l. and hence the values of b and Ti are well determined 
from AMscr{T). 

When ai = 0, the EPNJL model is reduced to the PNJL 
model. The results of the PNJL model are shown in Pig. |3]for 
comparison. The PNJL results cannot reproduce LQCD data 
particularly in the region T > 180 MeV. The slope of the solid 
and dot-dash lines in the region is thus sensitive to the value 
of ai- Namely, the value of ai is well determined from the 
slope. 

In Pig. |2] the solid and dot-dash lines are lower than the 
threshold mass 2Mj=i^n=o,a=o (dotted line). This guarantees 
that the and scr determined from the pole location 

in the complex-q plane agree with those from the exponential 
decay of (r) at large r. 

In the EPNJL model with the present parameter, the chi¬ 
ral susceptibility xii for light quarks has a peak at T = 163 
MeV, as shown later in Pig. fTTT a). This indicates Tg = 163 
MeV. The model result is consistent with LQCD data Tg = 
154±9MeVofRefs. Hi] . Por the deconhnement transition, 
meanwhile, the parameter Tq is adjusted to reproduce LQCD 
data on j'deconf^ already mentioned in Sec. |II] In fact, the 
Polyakov-Ioop susceptibility x<i>$ has a peak at T = 165 MeV 
in the EPNJL model, as shown in Fig. fTTT b). The model 
result X'^econf _ 265 MeV is consistent with LQCD data 


Fig. 3: Effects of T-dependent KMT interaction on pion and ao- 
meson screening masses. The solid (dot-dash) line denotes M.,r,scr 
(Mag,aci) calculated by the PNJL model with T-dependent coupling 
K{Tj. See Fig.[2]for LQCD data. 


J^deconf = 170 ± 7 MeV of Ref. jUl. 

Figure |4] shows T dependence of the renormalized chiral 
condensate defined by 




/,s — 








(39) 


and the Polyakov loop <1>. The present EPNJL model well 
reproduces LQCD data H for the magnitude of Z\; g in ad¬ 
dition to the value of Tg. The present model overestimates 
LQCD data for the magnitude of although it yields a re¬ 
sult consistent with LQCD for ydeconf overestimation in 
the magnitude of is a famous problem in the PNJL model. 
Actually, many PNJL calculations have this overestimation. 
This is considered to come from the fact that the definition 
of the Polyakov loop is different between LQCD and the 
PNJL model In LQCD the dehnition is I^lqcd = 

(trc Pexp[f Jp drA 4 (r, x)])/3, while in the PNJL model 
based on the Polyakov gauge and the mean held approxima¬ 
tion the dehnition is t^pNjL = ti'c exp[z(A 4 )/T]/ 3 , although 
both are order parameters of Z 3 symmetry dTII ; see for 
example Ref. as a trial to solve this problem. 

Now we investigate effects of T-dependent KMT interac¬ 
tion K{T) on M^ scr and Mag^sa- In Fig. |5] T-dependence 
of K{T) is switched off; namely, results of the EPNJL model 
with K{T) = K{Q) are shown. One can see that T- 
dependence of K{T) reduces the mass difference between 
and Mag,SCI signihcantly in a range 150 < T < 
180 MeV, comparing Pig. |5] with Pig. |2l At T = 176 
MeV where hrst-order chiral and deconhnement transitions 
take place, scr has a jump while Mag,sci has a cusp. Me¬ 
son screening mass is thus a good indicator for a hrst-order 
transition. 

In Pig. | 6 ] both the T dependence of K{T) and the entan¬ 
glement of Gs{'P) are switched off. Namely, the results of 
the standard PNJL model with a constant K are shown. The 
model cannot reproduce LQCD data, as expected. 

Figure [T] shows three types of EPNJL calculations for the 
mass difference AMsci{T). The mass difference plays a role 
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50 100 150 200 250 



T [MeV] 


T [MeV] 


Fig. 4: T dependence of Ai^s and The solid (dot-dash) line corre¬ 
sponds to results of the EPNJL model for Ai^s (^)- LQCD data for 
2-1-1 flavors are taken from Ref. 0. 



T [MeV] 


Fig. 5: Effects of T-dependent KMT interaction on pion and oq- 
meson screening masses. The solid (dot-dash) line denotes M,r,scr 
(Mag, bct) calculated by the EPNJL model with K{T) = iT(0). See 
Pig. Elfor LQCD data. 

of the order parameter of the effective restoration of U{1)a- 
The full-fledged EPNJL calculation (solid line) with both T- 
dependent K and the entanglement coupling Gg (^) well re¬ 
produces LQCD data, while the standard PNJL model (dotted 
line) with constant K largely overestimates the data. 

The present model has T dependence implicitly in Gs(^) 
through ^ and explicitly in K (T). As a model opposite to the 
present one, one may consider the case that K{T) = K(Q) 
and Gs has T dependence explicitly, i. e., Gg = Gg{T). We 
can determine GgiT) so as to reproduce LQCD data for 
however, this model overestimates LQCD data for AMga. 
Thus the present model is well designed. 


B. Meson susceptibilities 

The validity of K{T) is investigated by comparing LQCD 
data with the model results for meson susceptibilities 
(^ = 7r,^?7ns,(Tjis). LQCD data based on domain-wall 
fermions are available for two cases of pion mass M^r (0) 


Pig. 6: Effects of the entanglement coupling Gs(^) on pion and oq- 
meson screening masses. The solid (dot-dash) line denotes 
(Mag,sci) calculated by the standard PNJL model with constant K, 
i.e., the EPNJL model with K(T) = 7f(0) and ai = 0. See Pig. 
for LQCD data. 



T [MeV] 


Pig. 7: Mass difference AMsci{T) between pion and ao-meson 
screening masses. The solid, dot-dash and dotted lines denote re¬ 
sults of the EPNJL model, the EPNJL model with K{T) = K{0) 
and the standard PNJL model with K{T) = K{0), respectively. See 
Pig.[2]for LQCD data. 


at vacuum being the physical value 135 MeV and a slightly 
heavier value 200 MeV. In order to reproduce these values 
with the EPNJL model, we take mi = 5.68 MeV for the first 
case and 12.8 MeV for the second one. 

We consider the difference Att.oo = “ Xao^ or¬ 

der parameter of the effective [/(l)A-symmetry restoration. 
Ligure [8] shows T dependence of ATr^ag/T'^ for two cases of 
M^(0) = 135 and 200 MeV. Since the x™® have ultravio¬ 
let divergence, they are renormalized with the MS scheme in 
LQCD. Lor this reason, one cannot compare the LQCD data 
with the results of the EPNJL model directly. We then multi¬ 
ply the model results by a constant so as to reproduce LQCD 
data at T = 139 MeV for the case of = 135 MeV. The 

model results thus renormalized well reproduce LQCD data 
for any T in both cases of M7r(0) = 135 and 200 MeV. 

A similar analysis is made for T dependence of LItt.ct = 
xr" - xTZ “‘I = xZb - Aao" *at are related to 
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Fig. 8: T dependence of the difference ZiTr.uo between tt and ao 
meson susceptibilities for two cases of M,r(0) = 135 and 200 MeV. 
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Fig. 9: T dependence of and An,ao for M-„(Q) = 135 MeV. 


SUi^{2) X S'[/r( 2) symmetry. Figure l9l (fT0l i shows T depen¬ 
dence of and for M^(T) = 135 (200) 

MeV. In both the hgures, the EPNJL model well reproduces 
T dependence of LQCD results. The present model with the 
K{T) of (IHl is thus reasonable. 


C. The order of chiral transition near the physical point 

Finally we consider the order of chiral transition near the 
physical point = (6.2[MeV], 175[MeV]) in 

the mi-nis plane. First we vary mi from 9.9 to 0 MeV with 
TOs hxed at 175 MeV. 

Figure [TT] presents T dependence of the chiral suscep¬ 
tibility Xu for light quarks and the Polyakov-loop suscep¬ 
tibility x<z><f in three points, “simulation point (S-point)” 
of {mi,ms) = (9.9[MeV], 175[MeV]), “physical point 
(P-point)” of {mi,ms) = (6.2[MeV], 175[MeV]) and 

“light-quark chiral-limit point (C; point)” of {mi, mg) = 
(0[MeV], 175[MeV]). In general, Tc and J’^econt determined 
from peak positions of xu and x^^ depend on mi and mg. 
However, as shown in panel (a), the Tc thus determined is 
163 MeV at S-point and 160 MeV at P-point, and hence the 


Fig. 10: T dependence of zl,r,o- and Arj^ao for M^{0) = 200 MeV. 


value little varies between the two points. At C/ point, xii di¬ 
verges at r = Tc = 153 MeV. The chiral transition is thus 
second order at C/-point at least in the mean-field level. This 
result suggests that the effective t/(l)A restoration is not com¬ 
pleted at T = Tc. This suggestion is supported by LQCD 
data at S-point in Fig. [Tjwhere AMsci{Tc) is about a half of 
AMscr{0). 

As shown in panel (b), mi dependence of even 

smaller; namely, ydeconf _ ^55 MeV for S-point and C/- 
point and 163 MeV for P-point. At C/-point, has a sharp 

peak at T = 153 MeV. It is just a result of the propagation 
of divergence from xii to x^^ if^ . and never means that a 
second-order deconhnement takes place there. 

Next, both mi and mg are varied near P-point. Figure [T2| 
shows the value of log[xH(Tc)] near P-point in the mi-mg 
plane. The value is denoted by a change in hue. Three second- 
order chiral transitions (solid lines) meet at ~ 

(0,0.726TOP^y") = (0[MeV], 127[MeV]). This is a tricritical 
point (TCP) of chiral phase transition. 


IV. SUMMARY 

In summary, we incorporated the effective restoration of 
[/(1)a symmetry in the 2h-1 flavor EPNJL model by introduc¬ 
ing a T-dependent coupling strength K{T) to the KMT inter¬ 
action. The T dependence was well determined from state-of- 
the-art 2h-1 flavor LQCD data on pion and oo-meson screening 
masses. To derive the meson screening masses in the EPNJL 
model, we extended our previous prescription of Ref. for 
2 flavors to 2h-1 flavors. The strength K{T) thus obtained 
is suppressed in the vicinity of the pseudocritical temperature 
of chiral transition. As a future work, it is quite interesting to 
clarify how the present suppression is explained by instantons. 

In order to check the validity of K{T), we analyze 
TT, ao, rjns, (Tns-meson susceptibilities obtained by state-of-the- 
art LQCD simulations with domain-wall fermions ll2^ . The 
EPNJL model with the K{T) of (|6ll well reproduces T de¬ 
pendence of LQCD data. The present model building is thus 
reasonable. 

Using the EPNJL model with the present parameter set. 



A^/T(LQCD)-— 

^,«„^(LQCD) ..... 

A^yr^ EPNJL, KiT) - 

EPNJL, K(T) . 

/W„=135 MeV 
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between pion and ao-meson is effective to determine T de¬ 
pendence of the KMT interaction. The EPNJL model with 
the present parameter set is useful for estimating the order of 
chiral transition at the light-quark chiral-limit point and the lo¬ 
cation of the tricritical point, since it is hard to reach the chiral 
regime directly with LQCD. 

The present model consists of <?-dependent four-quark in¬ 
teractions and T-dependent six-quark interactions. Mean¬ 
while, the importance of eight-quark interactions was pointed 
out in Ref. 11^ . since it makes the thermodynamic potential 



m, [MeV] 


Fig. 11: T dependence of (a) chiral susceptibility xii (b) 
Polyakov-loop susceptibility at S-point, P-point and C;-point. 
Here xu and X<p 4 are dimensionless and their definition is the same 
as in the LQCD formulation. Calculations are done by the EP¬ 
NJL model with the present parameter set. The dotted, dot-dash 
and solid lines stand for the results at S-point, P-point and Cj- 
point, respectively. At Cj-point, xii is divided by 10 and diverges 
atT = T^ = 153 MeV 


we showed that, at least in the mean field level, the order 
of chiral transition is second order at the light-quark chiral- 
limit point of mi = 0 and mg = 175 MeV (the physi¬ 
cal value). This result indicates that there exists a tricriti¬ 
cal point near the light-quark chiral-limit point in the mi-mg 
plane. We then estimated the location of the tricritical point 
as {mums) Ri (0[MeV], 127[MeV]). 

In conclusion, we present a simple method for calculating 
meson screening masses in PNJL-like models. This allows us 
to compare model results with LQCD data on meson screen¬ 
ing masses. Meson screening masses are quite useful to de¬ 
termine model parameters. In particular, the mass difference 


Fig. 12: Order of chiral transition near physical point in the mi- 
rria plane. The value of log[xii(7c)] is shown by a change in hue. 
Simulation point, physical point, light-quark chiral-limit point and 
tricritical point are denoted by S, P, Ci and TCP. The solid lines stand 
for second-order chiral transitions. 

bounded from below. Furthermore, it is reported in Ref. 
that current-quark-mass dependence of quark-quark interac¬ 
tions is effective to reproduce meson pole masses with good 
accuracy. Therefore, further inclusion of these interactions is 
interesting as a future work. 
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